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We study directed percolation at the upper critical transverse dimension d — 4, where critical 
fluctuations induce logarithmic corrections to the leading (mean-field) behavior. Viewing directed 
percolation as a kinetic process, we address the following properties of directed percolation clusters: 
the mass (the number of active sites or particles), the radius of gyration and the survival probabil- 
ity. Using renormalized dynamical field theory, we determine the leading and the next to leading 
logarithmic corrections for these quantities. In addition, we calculate the logarithmic corrections to 
the equation of state that describes the stationary homogeneous particle density in the presence of 
a homogeneous particle source. 

PACS numbers: 05.40.-a, 64.60.Ak, 64.60.Ht 
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I. INTRODUCTION 

Directed percolation (DP) [jj is an anisotropic variant 
of the usual isotropic percolation (IP) Q in which an 
effect or activity can percolate only along a given pre- 
ferred (longitudinal) direction. DP is perhaps the sim- 
plest model leading to self-affme fractals. It has many po- 
tential applications, including fluid flow through porous 
media under gravity, hopping conductivity in a strong 
electric field , crack propagation £| , and the propaga- 
tion of surfaces at depinning transitions in one dimen- 
sion p|. Moreover, it is related to self-organized critical 
models Q. 

Often the longitudinal direction is viewed as time and 
DP is interpreted as a spreading process. In this dy- 
namic interpretation DP has become famous as the the 
generic universality class for phase transitions between 
an active and an absorbing inactive state The per- 

haps most intuitive spreading process belonging to the 
DP universality class is the so called simple epidemic 
process (SEP). In epidemic processes, individuals (also 
referred to as particles and for simplicity assumed to be 
located on the sites of a <i-dimensional lattice) are either 
susceptible, infected or immune. At time t an infected 
particle can randomly infect any of its susceptible neigh- 
bors with a certain activation rate. At t + 1 the newly 
infected particles are capable of infecting their suscepti- 
ble neighbors and so on. With a certain deactivation rate 
any infected particle may become immune. Depending on 
the difference r between the activation and the deactiva- 
tion rate, the process is endemic or epidemic. For r < t c 
the process dies out after a finite time. For r > t c the 
process spreads over the entire lattice and approaches a 
homogenous steady state. The point r = t c marks a 
non-equilibrium phase transition. There are 2 basic vari- 
ants of epidemic processes. In the SEP there is a chance 
that an immune particle becomes susceptible again. The 
spatio-temporal patterns generated by the SEP are DP 
clusters 0, S 13 • If immune particles remain immune at 



all times, one has the so called general epidemic process 
(GEP). The clusters of immune particles generated by 
the GEP are IP clusters. 

Since DP represents the generic universality class for 
phase transitions into absorbing inactive states, it fre- 
quently occurs that work on non-equilibrium phase tran- 
sitions addresses the question, whether a given system 
belongs to the DP universality class or not. Usually this 
is done via the fluctuation induced anomalous critical 
exponents. However, directly at the upper critical di- 
mension, i.e., in d — 4 transversal dimensions for DP, the 
leading scaling behavior is purely of mean-field-type and 
there are no anomalous critical exponents. However, fluc- 
tuations lead logarithmic corrections to the mean-field 
behavior. Just as the anomalous critical exponents, the 
logarithmic corrections can be used to decide if a given 
system belongs to the DP universality class. With the 
computer resources available today, numerical simulation 
on non-equilibrium systems explore more and more often 
high spatial dimensions. Simulations with reliable statis- 
tics of such systems in 4 dimensions are within reach 
today. Hence, we feel that it becomes important to know 
logarithmic corrections for the DP universality class. 

The leading logarithmic corrections are fairly easy to 
extract from the known renormalization group (RG) re- 
sults on DP. Astonishingly, this has not been done to 
date, at least to our knowledge. It has to be expected, 
though, that knowing the leading logarithmic corrections 
is not sufficient to obtain a decent agreement between 
theory and simulations. This expectation is based on 
the experiences that have been made for another system 
for which logarithmic correction have been studied in- 
tensively by numerical and analytical means, viz. linear 
polymers \Wi lllj . Numerical work on DP in d = 4 in 
progress seems to corroborate our expectation |12( . 

The aim of this paper is to derive analytically loga- 
rithmic correction for DP up to and including the next 
to leading order. We focus three dynamical observables 
that are well suited for investigation by numerical simu- 
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lations, namely the number N(t) of infected particles at 
time t generated by a seed at the spatio-temporal origin 
(x = 0, t = 0), their mean distance R(t) from the ori- 
gin (radius of gyration) as well as he survival probability 
P(t) of the corresponding cluster. Furthermore, we deter- 
mine logarithmic corrections for the DP equation of state 
(EQS) that relates the homogeneous particle density M 
of the stationary state to r and an auxiliary constant 
homogeneous particle source h. For dimensions below 4, 
this EQS is known to 2-loop order [l^. Its logarithmic 
correction have not been addressed hitherto. 

Of course, logarithmic corrections to dynamic quanti- 
ties like N(t) are not only relevant for DP; they are like- 
wise important for dynamic IP at the respective upper 
critical dimension 6. For logarithmic corrections to dy- 
namic IP we refer to Ref. [Tj]- Also, logarithmic correc- 
tion influence the static properties of IP clusters, like for 
example their various fractal dimensions, and their trans- 
port properties. For logarithmic corrections in static IP, 
see Ref. 

The outline of our paper is as follows. In Sec. ITT1 we 
briefly review the SEP as a dynamic model for DP. We 
sketch the renormalized field theory of the SEP and cite 
previous RG results. Furthermore, we explain by solv- 
ing the RG equation in d = 4 how logarithmic correc- 
tions arise in DP. In Section lTTTl we derive the logarithmic 
corrections for the aforementioned dynamic observables. 
Section IIVI treats logarithmic corrections of the mean- 
field equation of state. In Sec. we give a few conclud- 
ing remarks. Details of our diagrammatic perturbation 
calculation are relegated to an Appendix. 



II. A BRIEF REVIEW OF DIRECTED 
PERCOLATION AND ITS DYNAMICAL FIELD 
THEORY 

This section is intended to provide the reader with 
background on the dynamical field theory of DP and 
to establish notation. Moreover, it demonstrates how 
logarithmic correction emerge in the RG framework by 
solving the RG equation directly in d=4. 



A. Modelling directed percolation 

There are basically two complimentary approaches to 
model DP. The first approach is based on bond percola- 
tion and assigns a direction to the bonds. An example 
for this kind of model is the random resistor diode net- 
work, see, e.g., Refs. mm 13. 

In the other approach 
one models DP as a kinetic growth process fl9t I20L f2l| . 
viz. the SEP that we elaborated on in the introduction. 
Here, we will take the latter route. 

On mesoscopic scales it makes sense to describe the 
SEP in terms of the density n(x, t) of infected particles 
at time t and space coordinate x. It is well known that 
the Langevin equation (in the Ito sense) governing the 



time evolution of this density is given by |7J 

\~ x d t n{x, t) = V 2 n(x, t) - rn(x, t) - |n(x, i) 2 

+ C(x,i), (2.1a) 
C(x,i)C(x',t') = A _ V n(x, t)6(t - t')8{x - jsf) . (2.1b) 

The parameter r is essentially the rate difference men- 
tioned in the introduction and hence specifies the devi- 
ation from criticality. A represents a kinetic coefficient. 
C(x, t) is a Gaussian random field that subsumes reac- 
tion noise and otherwise neglected microscopic details. 
7TT stands for averaging over the distribution of the noise. 
The right hand side of Eq. I|2.1bjl goes to zero for van- 
ishing n(x, t) to enable the existence of the absorbing 
state. Contributions to Eqs. (|2.1() that are of higher or- 
der in the field or the derivatives turn out to be irrelevant 
in the sense of the RG. For example, a diffusional noise 
contribution can be neglected. 

Langevin equations are fairly intuitive and thus pro- 
vide a natural starting point in the mesoscopic descrip- 
tion of stochastic processes. Prevalent alternative forms 
of mesoscopic description are Fokker-Planck equations 
as well as dynamic functionals [23. I23I |24| . Dynamic 
functionals, also known as response functionals, are best 
suited for the application of field theory and RG ideas. 
This is the form of description that we will use here. The 
dynamic functional J for DP has been known for a long 
time 0, H] ■ After exploiting a rescaling form-invariance 
that allows to equate g and g' , J can be written as 

3 = J d d xdt\~s(\- l d t + {T- V 2 ) + |(s-S))s. (2.2) 

The order parameter field s(x, t) is proportional to par- 
ticle density n(x,t). s(x,t) is the response field corre- 
sponding to s(x, t). Time reflection, also known as dual- 
ity transformation, 

S(x,i) «— » -s(x,-t), (2.3) 

is a symmetry transformation of the dynamic func- 
tional 7]. This, however, is merely an asymptotic sym- 
metry that holds provided that irrelevant terms are ab- 
sent. When applying RG ideas to calculate leading scal- 
ing scaling properties or logarithmic corrections, one has 
to make sure that this symmetry is preserved. On the 
other hand, if one is interested in corrections to scal- 
ing stemming from irrelevant contributions to the func- 
tional 12.2fl . then one has to admit composite fields that 
break the symmetry H2.3(l . 



B. Renormalization and scaling 

The great virtue of response functional J is that it al- 
lows for systematic perturbation calculation in the cou- 
pling constant g that resembles many features of, and 
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allows to glean techniques from, the well established di- 
agrammatic perturbation treatments of equilibrium crit- 
ical phenomena. The most economic way to actually do 
these calculations is to use dimensional regularization in 
conjunction with minimal subtraction (minimal renor- 
malization). For background on these methods we refer 
to Refs. [2(j. An appropriate renormalization scheme 
is 

s->s=Z 1/2 s, S^S = Z 1/2 S, (2.4a) 
A^ \=Z- 1 Z X \, T^T = Z^ 1 Z r T, (2.4b) 
g -» g =Z- l Z- 1 ' 2 Z 1 J 2 G- 1 ' 2 ^/ 2 u 1 ' 2 , (2.4c) 

where the symbol ° indicates unrenormalized quantities. 
The factor [i e l 2 , where /i is an arbitrary external in- 
verse length scale and e = 4 — d measures the deviation 
from the upper critical dimension makes the renormal- 
ized coupling constant u dimensionless. The quantity 
G e = r(l + e/2)/(47r) d / 2 naturally appears in the com- 
putation of Feynman diagrams and is included here for 
later convenience. In minimal renormalization the crit- 
ical point value r = r c is formally set to zero by the 
perturbation expansion. In general, however, r c is a non- 
analytical function of the coupling constant. Thus, an 
implicit additive renormalization r — t c — > r is concealed 
in the minimal renormalization procedure. The renor- 
malization factors Z, Z\, Z T , and Z u are known to 2-loop 
order 0,0. 

The critical behavior of the Green's functions G„.n = 
([s]"[s]")( cum \ where (• -.)( cum ) denotes the cumulants 
with respect to the statistical weight exp(— J), is gov- 
erned by the Gell-Mann-Low renormalization group 
equation (RGE) 



Gn,h{{r,t};T,u;\,n) = 



with the RG differential operator T>„ given by 



/z<9 M + AC<9 A + Tnd T + j3d u . 



(2.5) 



(2.6) 



The Wilson functions featured in the RGE are known to 
2-loop order 0,0, 
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(2.7a) 
(2.7b) 
(2.7c) 



fi = -eu + ^ - (169 + 106 In t) ^ + O (u 4 ) . 

(2.7d) 

The RGE can be solved by the method of character- 
istics. The idea behind this method is to consider all 
the scaling parameters as a function of a single flow pa- 
rameter I. One sets sets up characteristic equations that 



describe how the scaling parameters transform under a 
change of /. The characteristic for the momentum scale 
fj, is particularly simple and has the solution jx{l) = fil, 
i.e., a change of I corresponds to a change of the exter- 
nal momentum scale. With help of the solution to the 
remaining characteristics one obtains 

G„, fi ({x, t}; r, u- A, = [^l) d Z{l)] (n+ " )/2 (2.8) 
x G„,f*({^x, Z x (l)(lfi) 2 Xt}; Z t {1)t/{vI) 2 ,u{1)\ 1, 1) • 

as a solution to the RGE. At this stage the scaling so- 
lution 12.8fl is still rather formal since Z(l), Z\(l), Z T (l) 
and u(l) require specification. Below the upper critical 
dimension, these quantities display power law behavior 
described by the well known critical exponents of the DP 
universality class. Directly in d = 4, their depend loga- 
rithmically on I and hence their behavior is qualitatively 
different from the lower dimensional case. 



C. General form of the logarithmic corrections 

Now we will state and solve the characteristics directly 
for d = 4. The Wilson functions cited in Eqs. I|2.7|l are 
central ingredients of the characteristics. To make the 
notion more economic, we will write /(it) = /o + f\u + 
+ ■ • • with / standing ambiguously for 7, £, k, and 
(3. We expect the meaning of the coefficients /o, /1 and 
so on to be evident. 

The characteristic for the dimensionless coupling con- 
stant u is given by 



, dw 
~dT 



l- 



(3{w) 



(2.9) 



where we abbreviated w = u(l). Solving this differential 
equation for e = yields 



/ = l(w) = IqW 



exp 



■J- + OH 

j3 2 w 



(2.10) 



where Iq is an integration constant. The remaining char- 
acteristics are all of the same structure, namely 



d\nQ(w) 



(2.11) 



Here, Q is a placeholder for Z, Z T , and Z\, respectively. 
q is a wildcard for 7, k, and £> respectively. Exploiting 
Id/dl = pd/dw we obtain the solution 



Q(w) 



] o w qi/02 exp 



(31 



w + 0{w 2 ) 



(2.12) 

where Qo symbolizes a non universal integration con- 
stant. 

The flow parameter introduced via the characteristics 
is arbitrary. This arbitrariness has an important virtue, 
viz. I can be chosen so that one of the relevant variables 



x, t, or t -1 effectively acquires a finite value in the scaling 
limit. In this paper we are interested in time-dependent 
quantities and hence we choose 



Z\(w)(ln) 2 \t = X , 



(2.13) 



where X$ is a constant of order unity. With this choice 
w and I tend to zero for A/z 2 t — > oo. Instead of using the 
original t, we find it convenient to use 



s = ^m(Vio)=~m(*/t ), 



(2.14) 



as our time variable. Here, to is a non universal time con- 
stant proportional to X . From Eq. (|2.10l) and Eq. I|2.12[) . 
specialized to Z\, we obtain 



s = w 1 — ao In w + 0(w) . 

The constant ao is given by 

foCi - 2ft 157 53 4 nQ7r „„ 

a = = In - = 0.976533 . 

2/3 2 192 96 3 

Finally, we find by using Eq. I|2.15(l 



w 



exp 



Ins 



a - 



O 



In s In s 1 



(2.15) 



(2.16) 



(2.17) 



for the dimensionlcss coupling constant as a function of 
time. 



III. LOGARITHMIC CORRECTIONS FOR THE 
NUMBER OF ACTIVE PARTICLES ETC. 

A. Number of active particles 

At criticality (r = 0), the number of active particles 
generated by a seed at the origin is related to the Green's 
function G\ i via 



N(t) 



d d x Gi i i(x, t] 0, it; A, fi) 



(3.1) 



Using the general scaling result (|2.8|l we can express the 
scaling behavior of N(t) as 

N(t) = Z{w) [ d d x(fil) d 



x G M (Z/ir, Z x {w)(ln) 2 \t; 0, w; 1, 1) 

= %)G 1 ,i(q = 0,lo;0,w;l,l) (3.2) 

with Z{w) given by Eq. H2.12|l specialized to Q = Z. 

Note that the Green's function in the last line of 
Eq. (|3.2|l depends on the renormalized coupling constant 
w. If we were interested only in the leading logarithmic 
correction to N(t) we could ignore this subtlety. The 
higher logarithmic corrections, however, are influenced 
by the specifics of the Green's function. For the second 
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FIG. 1: Self-energy E(q, t) at 1-loop order. 




FIG. 2: Dyson equation 13.411 to 1-loop order. 

logarithmic correction, we have to calculate Gi i i(x, t) to 
1-loop order. The diagrammatic elements required in this 
calculation are the Gaussian propagator 



G(q,t) =6(t) cxp[- X(T + q 2 )t] 



(3.3) 



where 9 stands for the step function, and the two 3-leg 
vertices \g and — Xg. In contrast to the calculation of 
for example critical exponents, it is not sufficient for our 
current purposes to consider Feynman diagrams with am- 
putated external legs because G\ t \ is determined by the 
Dyson equation 

Gi,i(q,i) = G(q,i) + f dt' f dt"G(<i,t-t') 
Jo Jo 

x E(q,t'-i")G(q, *") + ••• . (3.4) 

Here, E(q, t) stands for the self-energy that is given to 
1-loop order by 



E(q,t) = -^- J G(p,t)G(q-p,t) 



(3.5) 



where J stands for (2tt) d ^ 2 J d d p. Diagrammatic repre- 
sentations of the self-energy and the Dyson equation can 
be found in Fig. ^ and Fig. |3 respectively. For details 
of our calculation we refer to Appendix IA II Eventually 
we find 



Gi,i(q = 0, X ; 0, w; 1, 1) = 1 + A N (X )w + Q{w 2 ) . 



where An(Xq) is an amplitude given by 

1 



A N {X ) 



■\z + x) 



(3.6) 

(3.7) 

+ C'e 



Here, we used the shorthand notation Z = ln(2A"o) 
with Ce being Euler's constant. 

Now we know all the ingredients that contribute to the 
leading and the next to leading logarithmic correction to 
N(t). Collecting from Eqs. fFT^i and we find 



N(t) = No(w~ 



B 



B N ) 

,1/6 



1/6 



N 



1 - 



exp(— cnw + 0(w 2 )) 
6jv In s + cn 



O 



r In s In s 1 



(3.8) 
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where No is a non universal constant, Nq is a non univer- 
sal constant slightly different from N , and Bjy = 6Ajy. 
The first row of (|3.8|1 and the result (|2.15(1 constitute 
a parametric representation of the tuple (N, s) that is 
suitable for comparison to numerical simulations. The 
parametric representation has the advantage that it rep- 
resents a nicely systematic expansion in terms of the cou- 
pling constant w. The second row of (|3.8|) shows the 
more traditional form. The constants &at, cjv, and Bn 
are given by 



b N = 



7T 



157 



CjV 



1152 

72 

= 0.102113, 



R 71 



53 , 4 
In - 

576 3 

25 



1152 



0.162755. 

161 , 4 
m - 

576 3 



B N = -(Z 



1 



(3.9a) 

(3.9b) 
(3.9c) 



Equations (|2.14|) and (|3.8|l show directly that we may 
eliminate the arbitrary constant Z by a rescaling of the 
non-universal time constant to- However, we keep this 
constant in our formulas because it can mimic higher, 
neglected powers in our expansions, i.e., it represents a 
further constant that can be fitted to simulation results. 



B. Radius of gyration 

The mean square distance of the active particles from 
the origin, also known as their radius of gyration, is de- 
fined by 



R(t) 



_> /d d a;x 2 Gi,i(x,t) <91nG M (q, i) 



2d J d d xG!s(x,t) 



dq 2 



The general scaling solution (|2.8|l implies that 



q=0 

(3.10) 



dlnGi.ittZ/^q, Z x (w)(l/i) 2 \t; 0, w; 1, 1) 



Combining Eqs. I|3.11|l . (|3.12() and l|3.13|l as well as the 
solutions to the appropriate characteristics we obtain 



q=0 



_ 2 51nGi i i(q,X ;0,iu;l, 1) 



(3.11) 



q=0 



for t = 0. Equation (|3.11|) shows that the radius of gyra- 
tion is, like N(t), affected by the dependence of G^i on 
the renormalized coupling constant w. Here, however, we 
need the part of G\^\ that is quadratic in the momentum. 
In Appendix I A II we calculate that 



dq 



2 lnGi,i(q, X ;0, w; 1, 1) 



q=0 



= X (l + Ar{X q )w + 0(w 2 )) . 
The amplitude that appears here is given by 

A R (X ) = ±(Z-1), 



(3.12) 



(3.13) 



t^R 2 = Rl( W - l + B 
= R' 2 (s + B R 



1/12 , 

R ) expl-c R w 



0(w 2 )) 



1/12 



b R In s + c R 



O 



In s In s 1 



(3.14) 



with i?Q and R' 2 being non-universal amplitudes. 6r, cr, 
and Br = 12 Ar are constants that have the values 



cr 



ao 157 
12 ~ 2304 

6Jk 

P 2 02 2304 
0.0438136, 



53 4 

In - = 0.0813777, 

1152 3 



07 



59 , 4 
Tl52 ln 3 



Bn=" 7 (Z-l) 



(3.15a) 

(3.15b) 
(3.15c) 



C. Survival probability 

Recently it was demonstrated that the survival proba- 
bility P(t) of an active cluster emanating from a seed at 
the origin can be expressed in terms of the order param- 
eter and the response field as [23] 

P(t) = - km (c- kM S(-t)) , (3.16) 

k — >oo 

where J\f = J As(x,0). This formula is fundamental 
in that it relates the survival probability unambiguously 
to the fields inherent in the dynamic functional J . In 
actual calculations, however, the term exp(— kAf) has to 
be incorporated into the dynamic functional and one is 
led to 



Jk = J+ / dtk{t)M(t) 



(3.17) 



instead of the original J . Here, kit) = kS(t) is a source 
conjugate to the field s. Having introduced Jk 1 one can 
write 

P(t) = - lim (S(-t)) k 

k — >oo 

= — Go,i(— t, t, k = oo,u; X, fi) , (3.18) 

where (•••)* denotes averaging with respect to Jk- Note 
that the explicit term exp(— kAf) is gone. 

At this point we find it worthwhile to annotate an 
interesting implication of the time reflection symme- 
try (|2.3|) . Due to this symmetry the survival probability 
P(t) is identical to the mean particle density p(t) — (s(t)) 
of the dual process starting with a fully occupied initial 
state, p(0) — oo. 

To avoid tadpoles in our perturbation calculation, we 
carry out the shift s — > s + M so that (s) = is restored. 
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a) 



t 

-©- 



FIG. 3: (a) The new vertex -X 2 gM(t) and (b) the 1-loop 
tadpole diagram T(t). 



Go.i is then nothing but M. The entire procedure leads 
to the new response functional 



Jk 



j 1 . nil I a.;- U -1 ^ 

Ms 2 



(r - gM - V 2 



+ §(*-*))« + 



X 9 »~ r „ 2 



d t M + XrM - yM 2 



k)s 



(3.19) 



The diagrammatic elements implicit in comprise 
the two vertices encountered in Sec. IIII Al In addition, 
there is a third vertex, viz. — X 2 gM(t) as depicted in 
Fig-El 1 - The Gaussian propagator for the new functional 
has to be determined from the differential equation 



X^dt+r-gM^ + q 2 G{q,t,t') = X- x 6{t-t') . (3.20) 



To avoid tadpoles, M(t) must satisfy the differential 
equation 

d t M{t) ~ XrM(t) + yM (t) 2 - k(t) + f(t) = . (3.21) 

At 1-loop order, the tadpole T(t) is given by the diagram 
shown in Fig. [3t>- Upon solving the differential equa- 
tions i|3.20[l and i|3.21[) at the mean-field level we find 
that the modified Gaussian propagator reads (for details 
see Appendix I A 21) 



G (q,t,t')=&(t-t') 



Kpj-t) 
M-t>) 



exp 



X(T-q 2 )(t-t') 



where 



K (t) 



2t 



(e ATt - l) 



(3.22) 



(3.23) 



Having the modified Gaussian propagator at our disposal 
we are in the position to calculate the diagram depicted 
in Fig. |3Jd. This calculation leads eventually to 



G ,i(-X , 0, k = oo, w; 1, 1) 
cx uT 1/2 (l + A P (X )w + 0{w 2 )) 

with the amplitude Ap(Xq) reading 



A P (X ) = 1(Z- 3 - 



(3.24) 



(3.25) 



For details on this calculation we refer to Appendix I A 21 
Little further work is required to extract the logarith- 
mic corrections to the survival probability. Recalling the 
scaling form (|2.8(l and our choice for the flow parameter 
we deduce that, for r = 0, 

P{t) = -Z(w) 1 / 2 { f xl) 2 GoA-X , 0,oo, w; 1,1). (3.26) 

Inserting our results for the characteristics as well as 
Eq. I|3.24jl combined with l|3.25jl we obtain 

tP{t) = Patw- 1 + B P ) 1/2 exp(-c P w + 0(w 2 )) 



Pn(s + B f 



,1/2 



bp In s + cp 



O 



In s In s 1 



(3.27) 



Pq and Pq are non- universal factors. The constants bp, 
cp, and Bp — 2Ap are given by 



a 157 53 4 
<>,- = 4 = — + — In - = 0.488266 , 



cp 



B T 



2 384 

2C2 ~ 72 +/3. 
2p 2 ' ' 

0.00724268 , 

3 

2 



(3.28a) 



7i ~ 2Ci 



7 

384 



17 , 4 
m - 

192 3 



(3.28b) 
(3.28c) 



The constant Z might be eliminated by the same rescal- 
ing of the non-universal time scale to as discussed above. 



IV. LOGARITHMIC CORRECTIONS TO THE 
EQUATION OF STATE 

A. General considerations 

It is well known (see e.g. [2(J) that the generating func- 
tional W [ J, J] of the Green's functions 



G„, s ({x,t},{x,t}) = 



S n + h W\J, J] 



{8J(x,t)}"{8J(x,t)y 

(4.29) 

is related to the dynamic free energy functional T [s, s] 
by the Legendre transformation 



T+W = d d xdt[s(x,t)J(x,t) + S(x,t)J(x,t)] , (4.30) 



with 5W/SJ = s, 5W/SJ = S, 5T/5s = J and ST/SS = 
J. r is of great importance in diagrammatic perturba- 
tion theory because it is the generating functional of the 
irreducible vertex functions. The dependence of T upon 
the coupling constant g can be written in the form 



T[s,s;g] = g 2 $[gs,gs;u] 



(4.31) 
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The expansion of the functional &[(p,<p;u] into a series 
of u = G e /i~ £ g 2 yields the loop-expansion. The zeroth- 
order term g~ 2 $[<?s, gs; 0] is nothing else then the re- 
sponse functional J 12.21) itself. Hence, J constitutes 
the mean field part of the dynamic free energy. From the 
RGE (|2.5|) for the Green's functions it follows that the 
RGE for the renormalized dynamic free energy is given 
by 



d d xdt[s{^t) 



<5s(x, f) 



+ S(x, t)— - 



(5s(x, t) / 

x r[s, s; r, it; A, /x] =0. (4.32) 

Exploiting the findings of Sec. Ill CI the solution the 
RGE (|Q2l) is found to be 

r Is, s; r, u; A, /i] = (4.33) 

T [ZH^ 1 / 2 !, ^(u;)- 1 / 2 *; Z t («;)t, to; Z A (w)A, Z/x] . 

Now we revisit the Langevin equation (|2.1a|) . If the 
simple epidemic process is supplemented by an additional 
constant particle source h, Eq. I|2.1a|l is modified to 

A -1 9tra(x, t) = V 2 n(x, t) - rn(x, t) - |n(x, i) 2 

+ /i + C(x,t). (4.34) 

The extra term in the Langevin equation induces an extra 
term in the response functional (|2.2() . 



J^J h = J- J d d xdt A h s(x, i) . 



(4.35) 



Hence, the source J is shifted by Xh, and, to obtain the 
true generating functional of the irreducible vertex func- 
tions, we have to translate the field s by its mean value 
M. Then it is not difficult to see LUI that 



Xh = 



ST 



(5s(x, t) 



ri,o(Af,r) 



S=0.s=M 



'^ifiigM, t, u) 



(4.36) 



constitutes the EQS that relates the particle source to 
a constant mean particle density M for a given r. By 
simple dimensional considerations we find that <&io obeys 
the scaling relation 

$1,0 (gM, t, u; A, jtt) = A,u 4 $i,o (gM/n 2 , t/^ 2 ,u; 1, l) 

=:\ f i 4 F(gM/f J , 2 ,r/ l j 2 ,u). (4.37) 

A one-loop calculation ^| yields 

(x + y) 



F(x,y,u)/x = (y + x/2) 
+ 0(u 2 ) 



ln(x + y) — 1 

(4.38) 



Exploiting now Eqs. I|4.34ll and (|4. 37|t we get the result 

h = x /Gj^(l^) 2+d/2 Z x {w)Z(w)- 1/2 (4.39) 

x F(^w/G e Z(w)- 1/2 M/(l^ 2 , Z t (w)t/(1^) 2 ,w) , 

for the general scaling form of the EQS. Here, I = l(w) 
is given by Eqs. I|2.9|) and l|2.10|l . 



B. Behavior at the upper critical dimension d = 4 

In four dimensions we have G £= q = l/(47r) 2 . Following 
our work for d < 4 in U| we will cast the EQS in a 
parametric form. To this end we make the ansatz 



r/fi 2 = R(l - 



4ttM/h 2 = f M (w)R0 (4.40) 



so that R — corresponds to the critical point. The 
parameter 9 describes the crossover from the absorbing 
to the active phase. The source is zero for 9 = and 
9 = 9 > 1. We expect that 9q = 2. After inserting the 
ansatz (I4.40|) into the scaling form of the EQS (|4.39l) we 
choose the parameter w so that 



r 2 Z T (w)R = c, 



(4.41) 



where c is a convenient dimensionless positive constant. 
Defining 

p(w) = Z T (w)- 1 Z(wy 1/2 w 1/2 f M {w) (4.42) 

we arrive at 

4tt/i/^ 2 = R 2 Z x (w)Z T (w) 2 (wZ(w))- 1/2 

xc- 2 F(cp(w)9,c(l-9)) 1 w). (4.43) 

Next we determine the function p(w) so that 
F(cp(w)6, c(l — 9)), m) is analytic in 9 and has an ex- 
pansion in to. Using the 1-loop result l|4.38|l we obtain 
readily 9 = 2 + 0(w 2 ) and 

p(w) = 1 + w{l - lnc)/4 + 0(w 2 ) . (4.44) 

It follows that 

2c- 2 F(cp(w)9,c(l-9)),w) = 9(2~9) + 0{w 2 ) . (4.45) 

Using Eqs. (|2.10() . (|2.12() and l|4.41|l in conjunction with 
a rescaling of the arbitrary dimensionless variable i?, we 
get 



r := —~\nR = w 

4 



Oi In w + 0(w) . 



(4.46) 



The constant ai is given by 



«! = ^— = W2 + In - = 0.851533 . (4.47) 

Exploiting Eq. (|2.15() we obtain for the dimensionless 
coupling constant as a function of r the asymptotic ex- 
pression 



w = r exp 



lnr r^(^ n r mr 1 



(4.48) 



Collecting our results we get finally the equation of state 
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in parametric form 
t/tq = R(l - 
M/M = R6(w- 



y) 1/3 exp(-c M w 



(4.49a) 



R9(r + y) 



1/3 



b M In r + c M 



O 



In r lnr 1 



h/h = R 2 0{2 - 6) exp(-w/6 + 0(w 2 )) 



R 2 6{2 



1 - 



1 



O 



/lnr \ 



6r v r 

The constants 6a/, cjvf, and are given by 



6m = y = 

K 2 



133 
576 



53 4 

In - = 0.283844, 

288 3 



= 0.079712. 



71 

768 



17 

384 



(4.49b) 



(4.49c) 



(4.50a) 



(4.50b) 



Equations l|4.48|l and (|4.49b|) show directly that the ar- 
bitrary constant y can be eliminated by a rescaling of 
the nonuniversal constants To, Mq, and ho- However, we 
keep y in our formulas for the same reasons for which we 
kept the non-universal constant Z in Sec. II I II . 

To the order we are working here it is possible to elimi- 
nate the parameter R completely. Exploiting Eqs. i|4.46[l . 
(|4.49all and (|4.49c|) we can express R in terms of t and 
h as 

R = (1 + 0(w)) j (r/r ) 2 + h/ho . (4.51) 

Using Eq. i|4.51ll we can recast our results stated in 
Eqs. (j4~I%|) and {O^ as 



7(r/T ) 2 + h/h = w 4a ^ 3 



exp 



— +0(w) 



(4.52) 



and 

T/T 

M/M 

h/h 
(M/M ) 2 



■ y) 1/3 [i + c MW + o( w 2 )] , 

(4.53a) 
-yy V3 [l + (c M -1/6)w 

(4.53b) 



2-6> 



0(w 2 )} 



Finally, we recast our results in yet another form. As 
we will see shortly, this form highlights an interesting and 
probably important intricacy. If we take only the leading 
(1-loop) terms in Eqs. H4.52J) and l|4.53(l into account, our 
results boil down to the EQS 



M 





1/3 



(4.54) 



that is appealingly simple in stature. M' = (3/4) ^Afo 
is another non- universal constant. Now we switch to the 
scaled variables 



X = 



\/h/ho 



Y = 



M/Mq 



1/3 



(4.55) 



By incorporating the corrections coming from the 2-loop 
order we obtain 



Y = (Vl + X 2 - V 
+ 



1 



ln(l +X 2 ) 



1 

127 



Hh/ho) 

4ai lnr + 12c A f - 1 - 



1/3 



X 



In 2 



In this formula the variable r is given by 



r - - \Hh/h )\ 



1 



ln(l + X 2 



Hh/ho) 



(4.56) 



(4.57) 



Note that this asymptotic EQS l|4.56|l is not only a re- 
lation between the scaling variables (|4.55|l but that it 
also includes a dimensionful non- universal constant hg. 
Only in the limit |ln(/i//i )| — > oo we obtain the mean 
field equation of state although with the logarithmically 
corrected scaling variables l|4.55|l . This intricacy may in 
deed be relevant for the explanation of simulation re- 
sults El. 



V. CONCLUDING REMARKS 

In summary, we have investigated logarithmic correc- 
tions to scaling in DP by using renormalized dynamical 
field theory. We calculated the leading and the next to 
leading logarithmic correction for for the number N(t) of 
active sites at time t generated by a seed at the origin, 
the radius of gyration R(t) of the corresponding cluster 
as well as its survival probability P(t). Moreover we de- 
termined the logarithmic corrections to the mean-field 
equation of state that describes the dependence of the 
stationary particle density M (r, h) upon r and an auxil- 
iary external homogeneous source h. 

Our result involve 2 non-universal scales. The dynamic 
observables depend on the non-universal time-scale to 
and our asymptotic expansions are valid for times t ^ to- 
Note that to may serve as a measure of quality for mi- 
croscopic models of DP with respect to their suitability 
for numerical simulations. The smaller to is for a given 
model, the less computer time will be required in order 
to get good statistics on the critical behavior of DP. Our 
results for the EQS define the non-universal scale ho asso- 
ciated with the auxiliary source. The EQS depends on ho 
unless the limit |ln(/i//io)| — > oo is reached. Only in this 
limit one can expect logarithmically corrected mean-field 
behavior of the EQS. This should and probably has to 
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be taken into account when numerical data on the EQS 
is analyzed. The existence of the non-universal scales 
to and ho can be regarded as 2 examples for Coleman's 
concept of dimensional transmutation in naively scale- 
independent field theories. In other words, to and ho are 
akin to the hadronization scale of quantum chromody- 
namics. 

From the experience one has with other systems at 
their respective upper critical dimension, in particular 
linear polymers in d = 4, we expect that logarithmic cor- 
rections are of clear significance with respect to numerical 
simulations of DP in d — 4. This expectation is corrob- 
orated by recent simulations 01 • The aforementioned 
experience and the fact that we went up to the second 
logarithmic correction make us confident that our results 
compare will well with simulations, perhaps even quan- 
titatively. We hope that our analytical estimates trigger 
an increase effort to determine logarithmic correction for 
DP numerically with high accuracy. 
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APPENDIX A: CALCULATION OF GREEN'S 
FUNCTIONS 

In this Appendix we outline our 1-loop calculations of 
scaling functions belonging to the Green's functions Gi ; i 
and Go,i- These calculations provide us with the am- 
plitudes An(Xo), Ar(Xq), and Ap(Xq) that enter the 
second logarithmic corrections for the dynamic observ- 
ables. 



1. The Green's function Gi,i 

Here we provide some of the details of the calculation 
that leads from the Dyson equation lj3.4|l to our result for 
Gi.i stated in Eq. H^.fijl . First, we carry out the momen- 
tum integration in the self-energy l|3.5[) . This step leads 
to 



E(q,t) 



(A 3 ) 2 exp(-A(2r + g 2 /2)) 



(Al) 



2(87rAi)<V 2 

Next, we substitute (|A1(I into (|3.4|1 . After an integration 
we obtain 

G 1)1 (q,t) = G(q,t) 1 



i(2Xfi 2 t) 



e/2 



8T(l + e/2) 
x / dx (1 - x)x~ d l 2 exp(-aa;) 



(A2) 



Here, we introduced the shorthand notation 

a=(r-£)\t. (A3) 

The remaining integral is calculated in dimensional reg- 
ularization. A subsequent £-expansion yields 

dx (1 - x)x~~ d l 2 exp(~ax) 



(1 + a) 



1 + Ein(a) + 1 - exp(-a) + 0(e) , 

(A4) 



where the entire exponential integral is given by 



Ein(x) 



dy 



1 - exp(-y) 



y 



E 

fc=i 



k\k 



(A5) 



The next step is to remove the e poles by employing 

the renormalization scheme (|2.4|l . We let G\ : \ — > G\ \, 

A — > A, r — * f , and use the 1-loop results 

u „ u , . _ . 

Z x = 1 + — + ••• , (A6a) 



Z = 1 
Z T = 1 



(A6b) 



4e ' 8e 

u _ 2u 

— + •••, Z u = 1 + — + ••• . 

2e e 

Expressing the bare quantities in terms of these renor- 
malization factors and their renormalized counterparts 



Gij — ZG\.\ — ( 1 + ^J^M 



A = Z^ZxX 



8e 



Af = Z^Z t \t = [1 + — JAr, 



(A7a) 
(A7b) 
(A7c) 

and using the intermediate results in Eq. QA2fl combined 
with Eq. I|A4(1 we obtain the renormalized Green's func- 
tion 

G M (q, t) = G(q, t) 1 1 + | [(1 + a) ( ln(2A^ 2 t) + C E 

— Ein(a)J + a + cxp(— a) 

(A8) 

Note that the e poles are indeed removed by our renor- 
malization. 

Two results important for the logarithmic correction 
can be extracted from (| A8|) . Upon setting a = 0we find 

G M (q = 0,A^ 2 < = A ;r = 0,^1,1) 
= l + l(z + l)w, 



(A9) 

and hence the amplitude An(Xo) as stated in Eq. (|3.6|l . 

Moreover, we get 

-i 9 



Xn 1 — InGi,i(q, A/i t = A ; r = 0, w; 1, 1] 



q 2 =0 

(A10) 



which leads to our result for Ar(Xq) given in Eq. I|3.13|l . 
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2. The Green's function Go,i 

Here we provide selected details on our 1-loop calcula- 
tion of Gq i as required in i|3.18fl . We start by solving the 
differential equations i|3.2U[l and l|3.21[l . The initial and 
terminal conditions for the fields necessitate the ansatz 



M(t) = -9(-t)K(-ty 



(All) 



The type of the source term, k(t) — kd(t) with k — ► oo, 
demands the initial condition K(Q) = 0. With this in- 
formation, the differential equation i|3.21[l can be trans- 
formed without much effort into the integral equation 



it 



XtI 



dt'i 



-Art' 



K(t') 2 f(-t') + JL 



2t 
(A12) 

At mean field level, the solution to Eq. i|A12|) is given 
by Eq. l|3.23[l . Inserting the corresponding Mo(t) — 
— Kq(— i) _1 into the differential equation (|3.20|) we find 
the modified Gaussian propagator as stated in Eq. I|3.22ll • 
Now to the computation of the diagram depicted in 
Fig. 03. After some intermediate steps we obtain 



K{tfT{-t) = ^K (t)- 2 

K (t') 3 exp[2Ar(t - 1 1 )] 
X [8^A(< - <')] d/2 ' 



dt' 



(A13) 



The further evaluation of Eq. (|A13|) is fairly straightfor- 
ward for r = 0. After e expansion we find 



K(tff(-t) 



Xg 3 (2Xty/ 2 ( 6 



16(47r) d / 2 V e ' 2 
Insertion of this intermediate result into l|A12(l yields 



(A14) 



K(t) = 



gXt 



u{2\ t i 2 tf/ 2 / 3 9j 
4T(1 + e/2) V e + 4, 



(A15) 



Next, we renormalize. Indicating the consistency of our 
previous steps, the appropriate combination of renormal- 
ization factors Z^^Z^ 1 = 1 + 3u/ (4e) + • • ■ cancels the 
e pole in (|A15|> . The renormalized K (t) reads 



K(t) = f 



l-^(ln{\ft)+C E -\ 



(A16) 



Exploiting Go,i(— t) = K(t) 1 and X/j, 2 t = X as well as 
recalling the definition of 2 we finally obtain 



2irX G , 1 (-Xn 2 t = X ;t = 0, w; 1, 1) 



(A17) 
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